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We report an improved measurement of the partial-rate CP asymmetry in — » K^n^ decays. 
The analysis is based on a data sample of 85 million BB pairs collected at the T(4S) resonance with 

the Belle detector at the KEKB e+e" storage ring. We measure Acp(kK''') = 0.07 tom tom, 
where the first and second errors are statistical and systematic, respectively; the corresponding 90% 

confidence-level interval is -0.10 < Acp(k^'^^) < 0.22. 



PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 



In the Kobayashi-Maskawa (KM) model 1], CP violation arises from a complex phase in the quark- mixing matrix 
of the weak interaction. This idea is strongly supported by the observation of mixing-induced CP violation at the 
B-factories Direct CP violation (DCPV) is also expected in the KM scheme and has been observed in the K 
meson system However, DCPV has not yet been observed in the B meson system. 

Charmless hadronic B decays can provide opportunities to observe DCPV |5||, 0, 0, ■ Many of these decays 
include contributions from both h ^ u tree and b s penguin diagrams and the interference between these two 
processes can produce a partial-rate CP asymmetry: 

^ r(B ^ /) - r{B ^ /) 
'''' r{B^ f ) + r{B^ f) 

2|At||v4p| sin(5sin0 
~ \At\'^ + \Ap\^ + 2\At\\Ap\ cos 6 cos (j)' 

Here T{B — > /) denotes the partial width of either a B^ or 5+ meson decaying into a flavor-specific final state / and 
T{B /) represents that of the charge conjugate decay; At and Ap represent the tree and penguin amplitudes; and 
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5 and (p stand for the CP-conserving and CP- violating relative phases, respectively, between At and Ap. In order 
to have a sizable Acp, both phase differences have to be non-zero, i.e. S ^ and (j) ^ 0, and the tree and penguin 
amplitudes should be of comparable size {\At\ ^ \Ap\). 

The decay B± K° TT^ is almost a pure b ^ s penguin process and, thus, no sizable asymmetry is expected in the 
context of the Standard Model (SM) [aliot. However, our previously published result, based on an analysis of a 29 

fb~^ data sample, was Acp(k^'^'^) = 0.46 ± 0.15 ± 0.02 0. An asymmetry of this magnitude cannot be explained in 

the SM, even with the inclusion of interference of the basic penguin amplitude with a large — > (i^='=7r'')ticc K^tt^ 
re-scattering process \vi\ . and would be an indication of a new physics contribution in the penguin loop 12]. It is 
important to verify whether the central value persists with improved precision. 

In this paper, we report an updated measurement of the partial-rate CP asymmetry in K^tt^ decays based 

on a 78 fb~^ data sample collected at the T(4S) resonance, corresponding to 85.0 ± 0.5 million BB pairs, with the 
Belle detector at the KEKB e+e" storage ring This is approximately three times as much data as the 

sample that was used for the previous measurement and significantly improves the statistical precision. Throughout 

this paper, the partial-rate CP asymmetry AcpiK^T^^) is defined as 

A ('kK±^ N{Kl7:-)-N{Kl^+) 
■^^^(^ " ) = N{K%^-)+N{K%.+ y 

where NlK^n^) denotes the yield of B^ ^s'^^ decay and N{KgTT^) represents that of the charge conjugate mode. 

The Belle detector is a large-solid-angle spectrometer that consists of a three-layer silicon vertex detector (SVD), a 
50-layer central drift chamber (CDC), an array of threshold Cerenkov counters with silica aerogel radiators (ACC), 
time-of- flight scintillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) 
located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return located outside 
of the coil is instrumented to detect mesons and to identify muons (KLM). A detailed description of the Belle 
detector can be found elsewhere 13]. 

The analysis procedure is the same as described in Ref. _6j. Candidate B^ mesons are reconstructed using high 
momentum tt* and Kg mesons. For candidate tt* mesons, charged tracks are required to originate from the interaction 
region based on their impact parameters. 

In Belle, high momentum tt^ and mesons are distinguished by their associated Cerenkov light yield (-/Vp.o.) in 
the ACC and the ionization energy loss {dE/dx) in the CDC. These quantities are used to form a particle identification 
(FID) likelihood ratio TZ-^ ~ C.^/{C.^ + ^k), where £^ denotes the product of the individual likelihoods of A^p.c. and 
dE/dx for tt^ mesons; Ck is the product for mesons. For the TZt^ requirement used in this analysis, tt^ mesons 
are identified with an efficiency of 91% and there is a 10% misidentification rate. The efficiency and fake rate are 

estimated by comparing the'z?^ yields in a sample of £)*^-tagged*£'^ K^n^ decays before and after the appHcation 
of the high momentum FID requirements. A similar likelihood ratio that also includes the energy deposit in the ECL 
is used to identify electrons; positively identified electrons are rejected. Candidate mesons are reconstructed using 
pairs of oppositely charged tracks that have an invariant mass {ra.^.^) in the range 480 < to^tt < 516 MeV/c^. A 
candidate must have a displaced vertex and flight direction consistent with a K'g originating from the IF. 

Signal candidates are identified using the beam-energy constrained mass rabc = V -^bcam ~ P*b ^^'^ the energy 
difference A_E = E*g — -Ebcam, where i?bcam — 5.29 GeV and p*g and E*g arc the momentum and energy of the 
reconstructed B meson in the e^e~ center-of-mass frame. 

The dominant background comes from the e+e^ ^ QQ (q = u, d, s, c) continuum process; backgrounds from 
b ^ c decays are negligible because the momenta of the decay products are smaller than those of the signal Kg and 
TT^. We discriminate the signal from the qq background by the event topology. This is quantified by the Super-Fox- 
Wolfram (SEW) variable ja, 0] , which is formed from modified Fox- Wolfram moments that are combined using a 
Fisher discriminant into a single variable. The angle of the i3-meson's flight direction with respect to the beam 
axis {0b) provides additional discrimination. A likeHhood ratio TZs — Cs/{Cs + Cqq) is calculated, where Cs {C-qq) 
denotes the product of the individual SEW and Op likelihoods for signal {qq background). The probability density 
functions (FDFs) are derived from GEANT-based Monte Carlo (MC) simulations p/7| for the signal and rribc sideband 
(5.2 < mfcc < 5.26 GeV/c^) data for the qq background. We make a requirement on TZs that eliminates 88% of the qq 
background while retaining 73% of the signal. 

Signal yields are extracted from the Ai? distributions of events in the mbc signal region (5.271 < rribc < 
5.287 GeV/c^), separately for the K^tt^ and Kgir^ final states. The signal reconstruction efficiency ^3] is esti- 
mated to be 12% based on the MC. The AE distributions are fitted using a binned maximum likelihood fit with 
three components: the signal, qq background, and other charmless B decays, as shown in Fig. ^ The signal FDF 
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is modeled with a Gaussian distribution taken from the signal MC and calibrated using a —^ D^{^ K^Tr^)Tr^ 
sample where a similar reconstruction procedure is applied. For the qq background, the PDF is modeled with a 
second-order polynomial with a shape that is determined from the nibc sideband data. For other charmless B decays, 
the PDF is taken from a smoothed histogram of a large MC sample. (The enhancement in the lower AE region is due 
to charmless B decay modes with an additional tt meson.) Except for the signal peak positions, the same PDF shape 
parameters are used for both B^ and B^ samples. The signal peak positions are determined separately for the _B+ 
and B~ samples since a small systematic difference between the two samples is observed. (This is discussed below.) In 
the fit procedure, all of the PDF shape parameters are fixed and all the normalizations are free parameters. The signal 
yields are found to be N{K'^tt+) = 104.4 tH'l and N{K'^tt-) = 119.1 tlH, and the partial-rate CP asymmetry is 

determined to be Acp(K^Tr^) = 0.07 tom- 

The stability of Acp(k^tt^) as a fmiction of the the selection requirements is tested by varying the qq suppression 

requirement. As shown in Fig. 13 the value of AcpiK^T^^) is stable when this requirement is changed. 

Detector-based biases in Kgir^ reconstruction are investigated using a sample of inclusive, high momentum contin- 
uum — > KgTT^ decays, where the daughter particles are required to satisfy the same kinematic requirements and 
reconstruction criteria as used for the signal. Separate fits to the mass distributions, shown in Fig. |3| indicate 
that the signal mass resolutions for the 5+ and B~ samples are consistent, but there is a 1.0 ± 0.1 MeV/c^ difference 
in the mass peak positions. This difference in the peak positions is caused by a difference between the momentum 
measurement for high momentum negative and positive tracks that is attributed to a residual detector misalignment. 
After accounting for this difference in peak positions, Acp{D^ — > Kgir^) is determined and listed in Table |I| Here 

the sign convention in the definition of AcpiD"^ —>■ K^tt"^) follows that of AcpIk^t^"^)- The observed (2.0 ± 0.8)% 

asymmetry is treated as a possible bias, and —2.8% is assigned as a systematic error in the AcpiK^T^^) measurement. 

Possible biases in the B reconstruction are examined using a sample of B^ D^{-^ K^'k^)'k^ decays where 
the entire reconstruction procedure, except for the Kg reconstruction, is applied. Fits to the separate B'^ and B~ 
AE distributions, shown in Fig. ^ confirm that the resolutions are consistent between the two samples, while a 
3.2 ± 0.5 MeV difference in peak positions, due to the same effect that was found in the -^s""* sample, is 

observed. After accounting for the difference in AE peak positions, AcpiB^ D^n^) is determined and listed 
in Table The absence of an asymmetry indicates there is no bias. Biases in the high momentum PID and qq 
suppression are also examined by removing each of them in the Acp{B^ —^ D^tt^) measurement. The results are 
given in Table |l| No biases are observed. 

Possible asymmetries in the detector response and reconstruction for the qq background are checked using events 
in the nibc sideband region. The application of the entire reconstruction procedure confirms that the AE shapes of 
the _B+ and B^ samples are consistent, and no bias is observed, as indicated in Table The absence of TZ-^- and 
7?,s-related biases are confirmed in the same manner as for the B± ^'d^{^ K^ti^)-k^ sample. 

In order to study the sensitivity to the signal and qq background PDF shapes, each shape parameter is independently 
varied by its la error. In addition, the signal shape parameters are also estimated from the actual B^ ^s'^^ 
samples by allowing them to be free parameters in the fits. The uncertainty in the contribution from other charmless 
B decays is estimated from the change in the asymmetry by fitting the region of AE > —0.1 GeV without those 
decays. The resulting relative changes in asymmetries are added in quadrature giving the fitting systematics of -1-0.014 
and -0.006. 

Because of the difference between the results presented here and the sizable asymmetry in our previous measurement, 

the asymmetries of different data sub-samples are examined. Figure [S] shows AcpiK^T^^) for each data sub-sample 

together with AcpiD^ Kgir^) as a reference. The variation of AcpiK^T^^) is independent of that in AcpiD^ 
KgTT^) and is consistent with statistical fluctuations. 

The total systematic error in the Acp{K^t^^) is evaluated from the quadratic sum of the Kgir^ reconstruction 
bias and AE fitting systematics. Finally, the asymmetry 



A /irh^±\ _ n 07 +0-09 +0.01 
Acp[K TT ) - 0.07 _o,o8 _o,o3 



-0.08 -0.03 

is obtained and a 90% confidence level interval 

-0.10 < Acpi^K^TT"^) < 0.22, 



is set, where Gaussian statistics are assumed and the systematic error is added linearly. 
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In conclusion, we have measured the partial-rate CP asymmetry in K^tt^ with 85 million BB pairs collected 

on the T(4S) resonance at the Belle experiment. The resulting Acp{K^t^^) = 0.07 qs ^0 03 consistent with zero 

at the current level of statistical precision. The 90% confidence level interval — O.IG < AcpiK^Tr^) < 0.22 is set, which 
is consistent with other measurements This result has a statistical precision below 10% and supersedes our 

previous measurement jf|. We do not observe a significant partial-rate CP asymmetry in B^ K^n^ and attribute 

the sizable AcpiK^T^^) found previously in a much smaller data sample to a statistical fluctuation. 
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TABLE I; Summary of the detector-based bias tests. For tests other than those with the Kgn"^ sample, Acp values 
determined without the high momentum PID (TZtt) and qq suppression (TZs) requirements are also listed. 

Samples Acp [%] 

AT^TT* 2.0 ±0.8 

B± ^'i3^(^ A±7rT)7r± 0.6 ± 1.7 

w/o 0.0 ±1.5 

w/o 7^a 0.0 ±1.4 

^ A^TT^^ mtc sideband data 0.9 ± 1.3 

w/o 0.5 ±0.9 

w/oTZs 0.5 ±0.4 
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FIG. 1: The AE distributions for the — > K^iv^ candidates divided into B~ (left) and (right) samples. The fit results 
are shown as the solid, dashed and dotted curves for the total, signal and qq background, respectively; the hatched area indicates 
the contribution from other charmless B decays. 
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FIG. 2: Acp(K^n^) as a function of the signal efficiency of the qq suppression (TZs) selection. The horizontal line and hatched 

area indicate the Acp{K^''^^) value and its statistical error for the TZa requirement used in the actual measurement. Note that 
the statistical errors for the different data points are strongly correlated. 




FIG. 3: The mass spectra for the — > KgTv^ candidates separated into D~ (left) and D"*" (right) samples, where the 
kinematic requirements and daughter particle reconstruction are the same as used for the B^ Kg-K^ signal. The fit results 
are shown as the solid, dashed and dotted curves for the total, — > KgTv^ and combinatorial background, respectively. 
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FIG. 4: The AE distributions for B± ^ ■k^)tt candidates separately for the B (left) and (right) samples 

after application of the entire reconstruction procedure other than that for the Kg. The fit results are shown as the solid, 

dashed and dotted curves for the total, ^*L>\^ n^)TT and combinatorial background, respectively. The enhancement 
in the lower AE region contains backgrounds from B — > D*tt^ and D^p. 




Data sub-sample 



FIG. 5: Acp{K TV ) in each data sub-sample. The horizontal line and hatched area show the central value and the la statistical 

error of the Acp{K^''^^) result reported here. The solid points with the statistical error bars represent the Acp{K^''^^) result 
obtained for each data sub-sample; the open points show Acp{D^ — > K^n^). The sum of the three left-most points corresponds 
to the 29 fb~^ data sample used in our previous measurement. 



